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We describe the back action of microwave-photon detection via a Josephson photomultiplier 
(JPM), a superconducting qubit coupled strongly to a high-quality microwave cavity. The back 
action operator depends qualitatively on the duration of the measurement interval, resembling the 
regular photon annihilation operator at short interaction times and approaching a variant of the 
photon subtraction operator at long times. The optimal operating conditions of the JPM differ 
from those considered optimal for processing and storing of quantum information, in that a short 
T2 of the JPM suppresses the cavity dephasing incurred during measurement. Understanding this 
back action opens the possibility to perform multiple JPM measurements on the same state, hence 
performing efficient state tomography. 



I. INTRODUCTION 

Recently, many of the benchmark experiments of cav- 
ity quantum electrodynamics (QED) [1-4] have been re- 
produced with superconducting circuits [5-12] , which op- 
erate in the quantum regime via exchange of microwave- 
frequency excitations [13-17]. On these circuits, non- 
linear devices couple to the microwave-frequency modes 
of transmission lines via ordinary circuit devices such as 
capacitors or inductors [18-20], much as atoms couple 
to modes of a resonant electromagnetic cavity. A fixed 
number of these artificial atoms can be fabricated in a 
given circuit, and their energy levels and interactions are 
tunable both at fabrication and during the course of an 
experiment. For these reasons, circuit-QED (cQED) re- 
ceives attention both as a possible platform for scalable, 
universal quantum computing [21, 22] and for its ability 
to operate in regimes inaccessible by atomic cavity-QED 
[23-30]. 

While many of the tools available in cavity-QED are 
straightforward to reproduce in the circuit analogue, the 
detection of single microwave-frequency photons proves 
challenging. Traditionally, the lower cutoff frequency of 
photon counters is determined by the work function or 
band gap of a certain material, which is at a minimum 
in the infrared range for stable materials. There are cur- 
rently a few theoretical proposals for the construction 
of microwave photon counters [31-33], and recently it 
was demonstrated experimentally that a current-biased 
Josephson junction [34-39] can be used to count mi- 
crowave photons [40]. We refer to such a device as a 
Josephson photomultiplier (JPM), distinguishing it from 
a phase qubit [41, 42] because the optimal operating con- 
ditions for photon detection are different than those re- 
quired for storage of quantum information. 

Photon counters should be contrasted with amplifiers. 
While the former are sensitive to the intensity of the in- 
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coming radiation but not to the phase, the latter am- 
plify the quadratures of the signal. Even though com- 
mercial microwave amplifiers operate far from the quan- 
tum limit, researchers have recently demonstrated single- 
photon sensitivity in phase-preserving microwave ampli- 
fication [43-51]. By contrast, phase-insensitive photon 
counters have proven useful in quantum optics for recon- 
struction of the quantum state of light as, for example, 
in homodyne tomography. Workarounds to microwave 
photon counting based on linear amplification have been 
formulated [52, 53] and demonstrated [54-56]. 

As microwave photon counters, JPMs have an impor- 
tant application in efficient quantum state tomography 
of microwave photon states. Given that measurement 
by a JPM provides only limited information - a click 
indicates the presence of one or more photons - the 
post-measurement state still contains coherent informa- 
tion about the initial state. Following the idea of quan- 
tum regression this post-measurement state is connected 
to the pre-measurement state by back action operators. 
Hence, if the back action operators are known, repeated 
measurements on a chain of post-measurement states can 
provide additional information about the original state. 

In this paper, we theoretically model the back action 
of the proposed JPM, obtaining the precise relation be- 
tween pre- and post-measurement states. This knowledge 
may allow efficient state tomography [57, 58] including, 
for example, adaptive techniques [59], or any other ap- 
plication requiring knowledge of the post-measurement 
state. We note that our results are based on a very ab- 
stract model and thus extend to other detection schemes 
whereby a quantum two-level system strongly couples to 
a resonant linear oscillator, so long as the observable de- 
tection event involves incoherent tunneling from an en- 
ergy level of the two-level system (and not the oscilla- 
tor). For example, this situation applies to some setups 
in atomic cavity QED. We include realistic estimations 
of the energy dissipation and dephasing rates of a JPM, 
showing that operating the JPM in the regime of fast 
dephasing (short T 2 ) reduces the amount of dephasing 
incurred during measurement, and is thus advantageous. 

In the following section, we discuss our model for the 
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JPM. In section III, we discuss the formalism of process 
tomography used to characterize the back action of the 
JPM. In section IV, we give the back action both numer- 
ically and analytically in a variety of instructive and / or 
experimentally relevant regimes. In section V, we discuss 
how to extract the operating regime of the JPM by sim- 
ple tests with coherent light. Finally, in section VI, we 
discuss the optimal working conditions of the JPM for 
the purpose of cavity state reconstruction. 



II. PHYSICAL MODEL 

A JPM consists of a current-biased Josephson junction 
(CBJJ) [34-39] capacitively coupled to the microwave 
cavity of interest. The potential energy of a JPM , shown 
in figure 1, is 
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where I c is the critical current of the junction externally 
biased by current is the superconducting phase dif- 
ference across the junction, and 5>o = ^ is the magnetic 
flux quantum. We consider a junction biased such that 
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FIG. 1: This figure is a diagrammatic representation of the 
potential energy of a JPM as a function of the superconduct- 
ing phase difference, and shows the interaction with a mi- 
crowave cavity. 

the potential well contains only a few meta-stable states. 
All of these states can tunnel incoherently out of the po- 
tential well, but due to the exponential relationship be- 
tween tunnelling rate and barrier height, the tunnelling 
rate for higher energy states is several orders of magni- 
tude higher than that for the ground state. Note that, 
in analogy to the phase qubit, the current source can 
be replaced by a large, flux-biased superconducting loop 
[42,60]. 

Photon detection relies on an incident microwave pho- 
ton to transition the JPM to its first excited state. This 
transition is enhanced by pulsing the bias current, bring- 
ing the energy level splitting of the ground and first ex- 
cited state of the JPM on- or near-resonant with the mi- 
crowave cavity. Once the JPM reaches its excited state, 



it tunnels more rapidly out of the metastable state. This 
tunneling process is incoherent, resulting in a measur- 
able voltage pulse in the circuit that is interpreted as the 
detection of a single photon [36-40]. By comparison, a 
related method of cavity state reconstruction determines 
the number of cavity photons present by the frequency of 
coherent oscillation between the cavity and a phase qubit 
[24]. 

We assume that the JPM and cavity states are initially 
separable with the JPM in the lowest energy metastable 
state. Physically, this means that the JPM and mi- 
crowave cavity must be brought on resonance adiabat- 
ically with respect to the JPM's internal evolution, but 
non-adiabatically with respect to the cavity- JPM in- 
teraction. We describe the full system's Hilbert space 
with tensor products of single mode cavity eigenstates 
and three detector states, {|0) d , |l) d , |m) d }. The states 
|0) d and |l) d correspond to the ground and first excited 
metastable states of the JPM, while the measured state 
|m) d is an amalgamation of the many possible states that 
the JPM can tunnel into incoherently, producing an ob- 
servable output voltage. 

The coherent interaction between the cavity and the 
JPM as well as the relevant incoherent processes (for ex- 
ample, tunneling into the measured state, dephasing, and 
relaxation) are described by the quantum master equa- 
tion 

i(t) = sm\ 

= -ziH,m+Y,( J ^ J t - ^4 £(*)}) (2) 

where is the cavity- JPM system's density matrix. 
Here H is the Jaynes-Cummings interaction 
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where a and are the lowering and raising operator as- 
sociated with the cavity mode, a are the lowering and 
raising operators between the states |0) d and |l) d , and 
g is the coupling strength between the cavity and the 
JPM. Note that this Hamiltonian conserves total excita- 
tion number and does not couple photons coherently to 
the measured state, which simplifies the following analy- 
sis. 

A set of Linblad operators { J^} describe the relavent 
incoherent processes. 



J x = V=n(Hc®|m)<l| d ) 



(4) 



describes incoherent tunneling out of the excited 
metastable state leaving one less excitation in the cavity, 
and thus corresponding to photon measurement. Tun- 
neling out of the metastable ground state is described 

by 



Jo = V7b(Ic ® |m) (0| d ) 



(5) 

where 70 is the effective dark count rate since a mea- 
surement signal is produced without changing the num- 
ber of excitations in the cavity. We also include Linblad 
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operators to describe pure dephasing of the JPM over 
characteristic time T2 
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and energy relaxation from the excited state to the 
ground state of the JPM over characteristic time T\ 



J 3 = 
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(Hc®|o><i| d ). 
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In general the cavity decoheres as well [61, 62], but this 
happens slowly compared to other relevant timescales. 



III. PROCESS TOMOGRAPHY 

For tomography of the process of cavity state measure- 
ment by a JPM, we calculate the Liouville supermatrix 
T(t) generated by 

S = -i (h ® i - 1 ® ff) 

+ ( h ® h - ® 1 - ^ ® Jlh 



which transforms an arbitrary, vectorized initial cavity- 
JPM state to the solution of (2): 



i(t) = T(t)£(0) = e St £(0) 



(8) 



We then transform the Liouville supermatrix into the 
more commonly used ^-matrix representation [63, 64], 
which, for a given basis {E^}^ = q 1 of operator space 
satisfies 
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In the standard basis, E^i a ^\ = \a) where /j,(a,(3) = 

(N xa) + f3 and {la)}^^ 1 ^ s an eigenbasis of the full sys- 
tem's noninteracting Hamiltonian, the ^-matrix elements 
are simply 
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and can be obtained by a permutation of the Liouville su- 
permatrix elements: Xafi-yS — T a7 pg. In all calculations, 
we assume an initial state of the form 



m = p c (o) ® |o)<o| d) 



(ii) 



a product state of the initial cavity state, p c (0), and the 
lowest energy metastable state of the JPM. Preparation 
of this factorized state has already been described earlier 
(section II). 



We are interested in the back action of the JPM onto 
the cavity state conditioned on measurement outcome 
se {0,1, m}, 



P s c (t) = 



P'{t) 



(12) 



where P s (t) normalizes the cavity state by the probabil- 
ity of obtaining JPM final state \s),. Only incoherent 
tunneling into the measurement state is allowed by our 
model, so no coherent superposition between the mea- 
sured and non-measured state is possible. Therefore, 
ps=m gj ves cav ity state after detection of a photon, 
and in the case that no photon is detected, the cavity will 
be in a mixture of the states with s = 0,1. Each outcome 
is completely described by an off-diagonal c?^ av x c^ av 
block of the full x- ma trix, which is by itself a valid \~ 
matrix of the isolated cavity. We label these reduced 
X-matrices \ s , which completely describe the evolution 
of an arbitrary initial state p c (0), in the case of measure- 
ment outcome s. 



IV. ANALYTICAL AND NUMERICAL 
SOLUTIONS FOR THE X -MATRIX 

A. No Tunneling Model 

To understand the back action of photon detection, we 
first consider a simpler model where there are no incoher- 
ent processes, and measurement of a photon corresponds 
to projecting the JPM onto its metastable excited state 
|l) d (rather than |m) d ). In this model, eq. (2) reduces 

to = —i[H, £,{t)], and the cavity conditioned on mea- 
surement outcome s at time t can be expressed in terms 
of a single time-dependent back action operator B s {t) 
acting only on the Hilbert space of the cavity: 



B s (t)£(0)B s (ty 



(13) 



Furthermore, in this model, the B s (t) are straightforward 
to calculate explicitly: 



B\t) = (mm d 

00 

= -i E sm(gts/n) \n - 1) (n|, 



(14) 



B\t) ee <0|^)|0) d = Ecos(^)|n)(n| c . (15) 

The cavity is initialized in a superposition of n-photon 
Fock states, and when only coherent cavity- JPM interac- 
tion is included, each n-photon Fock state in superposi- 
tion will exchange a single excitation with the JPM at a 
Rabi frequency g^fn. Measurement projects the detector 
onto |0) d or |l) d , modifying the cavity with back action 
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B°(t) or respectively. From these operators we 

obtain the average detection probability 

P\t) = (B l \t)B\t)) Q 

= J2 P n(t)(n\Pc(0)\n) (16) 

n 

where P^(i) = sin 2 (gty/n) is the detection probability 
when the n-photon Fock state is initially prepared. 

It is possible, by averaging over repeated measure- 
ments at increasing time intervals, to distinguish Fock 
states and incoherent mixtures of Fock states by Fourier 
transforming the average detection probability -P 1 (t), as 
was demonstrated in ref. [65]; however, more sophisti- 
cated state tomography is required for resolving super- 
positions of Fock states. One approach is to displace the 
cavity state and reconstruct a convenient phase space de- 
scription of its initial state, repeating measurements and 
averaging at every point in phase space that is resolved 
[24]. Here we look for a quantum description of the mea- 
sured cavity state so that repeated measurements on a 
single input state can be used for a more efficient state 
tomography. 

From the behavior of (14), 

B 1 = {-igt)a + 0{yfagtf, (17) 

we can see that at short-times (and finite n), the mea- 
surement back action is proportional to the photon low- 
ering operator. Furthermore, at short times, different 
Fock states can be distinguished by their tunneling rate 
into the detector, T n = g 2 n. However, at longer times 
t > t n = 1/g^/n, the oscillations from different Fock 
states become out of phase and difficult to distinguish. 
This effect is described by the correction to the a opera- 
tor in (17). 

We note that back action is well-described by the low- 
ering operator when the interaction times t n (for the 
largest occupied n) are the largest time scales of the 
system. This happens, for example, in the free photon 
regime where coupling strengths are very weak [66] . Also, 
we will later show that in the richer model where detec- 
tion corresponds to tunneling out of the (l| d at the rate 
7i , the lowering operator is a good approximation to the 
back action for times shorter than the excited state tun- 
neling time 7f 1 . To detect low-energy microwave pho- 
tons, however, we need measurement times long com- 
pared to t n and 7^ 1 , and corrections to the lowering op- 
erator become important to fully understand the back 
action associated with photon detection. 

B. Full (Tunneling) Model 

In the full model, a detection event corresponds to 
the detector incoherently tunneling to the 'measurement ' 
state |m) d , which corresponds to the JPM observably 
tunneling from a metastable state to a near continuum 



of levels. The JPM does not reset to its initial state on a 
time scale comparable to the measurement interval and 
therefore can not resolve the total number of photons 
present in the cavity, only whether there is at least one 
present. Therefore, we do not expect the measurement 
back action on the cavity to be exactly the photon anni- 
hilation operator, but rather an operator of the form 

oo 

B m = ^2\n-l)(n\. (18) 

n=l 

which we refer to as the subtraction operator. We note 
that this back action can be used to separate the number- 
dependent part of the annihilation operator, 

a = BkN 1 ' 2 , (19) 

and was thus considered as the exponential of a quan- 
tum phase operator by Susskind and Glogower [67], but 
a currently more accepted unitary version was proposed 
by Pegg and Barnett [68]. 

While in the no-tunelling model the measurement 
back action reflects undamped Rabi oscillations between 
the cavity and JPM, we expect these oscillations to be 
damped by incoherent tunneling out of the metastable 
states of the JPM so that, when averaged over an entire 
measurement interval, the back action has the form of 
(18). Transitions occur from |n) to \n — l) photon Fock 
states, with no preference on the number n of photons 
originally present. However, because the initial tunnel- 
ing rate depends on the number of photons present, and 
it takes time for this averaging effect to occur, we expect 
that for measurement intervals short compared to t n and 
7-j -1 the back action will more closely resemble that of the 
usual photon-number resolving annihilation operator. 

To understand the distinguishing signatures of a and 
B m in the framework of process tomography, it is instruc- 
tive to examine the x 1 matrices corresponding to each. 
Both will have 

Xj-ijk-ik = Pjk (20) 

V j,ke{i,...,N-i} 

non-zero, corresponding to superpositions of \j) and 
photons transitioning to \j — 1) and \k — 1) photons. For 
a good photon detector, the number of excitation in a 
given Fock state is decreased by exactly one, therefore 
all other elements of x 1 are zero. When the back action 
operator is a, 

Pjk = VJk, (21) 

while for B m , 

frk = 1 (22) 

for all values of j,k e {l,...,N — 1}. In the following 
section, we numerically study the time-dependence of the 
/3jk in our full physical model, using the values of j3jk 
for known examples of back action models as a point of 
reference. 
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C. Numerical Simulations for the x 1 Matrix 

1. Bare JPM 

Here we present the ^-matrix elements numerically 
generated using the Liouville supermatrix approach, first 
in the case of a bare detector experiencing no dark counts, 
dephasing, or energy dissipation (70 = and T\ = T2 — 
00). In this case, the x^-matrix has the same nonzero 
elements as those for a and B m , labelled above as [3jk- 
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(a) Bare JPM Diagonal Matrix Elements 




(b) Bare JPM Off-Diagonal Matrix Elements 



FIG. 2: (a) Diagonal and (b) off-diagonal ^-matrix elements 



for a bare JPM as a function of time, where atj 
denned in equation 20. 
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The are plotted as a function of total measurement 
time (i TO ) in figure 4 for j, k £ {1, 2, 3, 4}. As is demon- 
strated clearly by the diagonal ctj = f3jj plotted in figure 
2(a), the ^-matrix elements show oscillatory behaviour 
at n-dependent frequencies, as in equation (14), g\/n. 
Similarly, the off-diagonal elements also show oscillatory 
behaviour, as can be seen in figure 2(b). 

In the long time limit, the diagonal elements all tend 
to unity as expected for a back action resembling the sub- 
traction operator, however, the off-diagonal elements do 
not. This additional dephasing can be explained by the 
uncertainty in time of the switching event, which is of the 
order 7^ . As the phase of the off-diagonal matrix ele- 
ments precesses with frequencies proportional to g 7 this 



uncertainty gives a spread of the phase of size 5/71. We 
will later see how decoherence reduces this uncertainty, 
thus reducing the amount of dephasing incurred by mea- 
surement . 



2. Pure Dephasing 

We now consider a JPM that experiences pure dephas- 
ing between its ground and excited states, as would be 
described by a master equation including the Lindblad 
operator J 2 of equation (6). In this case, the x 1 - ma tri x 
has the same non-zero elements as that for the bare detec- 
tor since the selection rules imposed by the conservation 
of excitation number are still valid. The bare detector 
and pure dephasing x^-ni&trix elements are compared in 
figure 3, where the dephasing time has been chosen such 
that 7jr = 10 71, deep in the strong dephasing regime. 
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(a) Bare JPM and Pure Dephasing Diagonal Matrix Elements 
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(b) Bare JPM and Pure Dephasing Off-Diagonal Matrix Elements 

FIG. 3: Comparison of the (a) diagonal and (b) off-diagonal 
X 1 -matrix elements of a bare JPM with one experiencing ad- 
ditional pure dephasing I/T2 = IO71 (marked with circles). 
In both plots, each curve's color indicates its row in the \ l ~ 
matrix, and in (b), distance from the diagonal is indicated by 
line style. 

As can be seen in figure 3, decreasing the value of T2 
has multiple effects. On the one hand, the photon trans- 
fer from the cavity into the detector is slowed down by 
decreasing T2, so at short time the ^-matrix elements 
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are smaller. Also, the coherent oscillations are damped 
as T2 has the effect of turning the coherent tunnelling 
between the cavity and the JPM into an incoherent pro- 
cess, similar to the crossover from strong coupling cavity 
QED to the Purcell regime [3]. In fact, a phase Purcell 
effect has been discussed in [69, 70]. This affects both 
the diagonal and off-diagonal % -matrix elements. 

Once the photon transfer efficiency is no longer a limit- 
ing factor, the asymptotic limit of the diagonal ^-matrix 
elements is not affected by T 2 . The off-diagonal elements 
saturate to a value that is set by measurement-induced 
dephasing which is lowered by short T 2 (as seen in figure 
4(a)). This reduction in measurement-induced dephasing 
is due to the fact that T 2 turns the coherent tunnelling 
between cavity and JPM into an incoherent process, and 
thus reduces the phase precession of the off-diagonal ele- 
ments and with it the uncertainty of these phases at the 
moment of measurement. Although the amount of de- 
phasing decreases with decreasing T 2 , the total measure- 
ment time required for the ^-matrix elements to reach 
their asymptotic value increases, as shown in figure 4(b). 



herent state input and a one-photon Fock state input. In 
figure 5, we see that for both input states, dephasing sup- 
presses the oscillations in detection probability exhibited 
by the bare JPM. These oscillations result from coher- 
ent excitation swapping (Rabi-type oscillations) between 
the cavity and JPM, and superpositions of an excitation 
in the JPM and in the resonator are subject to dephas- 
ing processes in the JPM. This pure dephasing turns co- 
herent JPM-cavity oscillations into incoherent resonant 
tunneling. 

In the long time limit, both the dephased and bare 
JPM detect a photon with the same probability. Thus, 
it is not necessary to aim at long T 2 values for a JPM as 
one would for a phase qubit [41, 42]. On the contrary, 
we see that the dephasing incurred by measurement is 
smaller at short T 2 , rendering it advantageous. A more 
detailed discussion about T 2 as an engineering parameter 
will be given in the end of the paper. 
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(a) Off-Diagonal Asymptotic Limits 




(a) Coherent State Detection Probability 




9 T ? 




(b) Off-Diagonal Asymptotic Time Scale 



(b) One Photon Fock State Detection Probability 



FIG. 4: Figure (a) shows the asymptotic limit of the off- 
diagonal x matrix elements as a function of T2, and figure 
(b) shows the timescale over which these asymptotic limits 
are reached. Figure (b) is not monotonic due to the coherent 
oscillations present for long T2. 

To further illuminate the effect of pure dephasing on 
a JPM we look at the probability of detection for a co- 



FIG. 5: Figure (a) shows the detection probability for a bare 
JPM and a JPM experiencing pure dephasing for a coherent 
state input (as described in section (VA) with \a\ — 0.5). 
Figure (b) shows the detection probability for a bare JPM 
and a JPM experiencing pure dephasing for a one photon 
Fock state input. 
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3. Energy Relaxation of the JPM 

We consider a JPM experiencing energy dissipation, as 
described by a master equation including the Lindblad 
operator J3 of equation (7) . In this situation, the JPM 
will also experience associated dephasing on a timescale 
T2 = 2Ti. The ^-matrix of a dissipating JPM has ad- 
ditional non-zero elements (in addition to the j3jk) at- 
tributed to a change of total excitation number. As 
the JPM can lose photons into an external heat bath, it 
is possible that multiple photons from the cavity might 
excite the JPM before a detection event occurs. The 
nonzero x 1 -matrix elements are of the form 

Pjk = Xj-lj+rk-lk+n (23) 

corresponding to the loss of integer < r < minjfc, j} 
photons before detection, fiy^ — j3jk from the previous 

(r) 

sections, and the diagonal elements are relabeled ay = 

In figure 6 we compare the a^p for an energy relaxation 
timescale of jr = 71 to the otj of a detector with infinite 
T\. Unlike the case of pure dephasing, energy loss from 
the JPM reduces the asymptotic value of these ^-matrix 
elements. We also note that the off-diagonal evolve 
exactly like those of a JPM experiencing pure dephasing, 
as shown in figure 3(b), but with an effective T 2 of 2Ti. 

The diagonal ^-matrix elements a*p are shown in figure 
7 for different values of r. 




FIG. 6: This figure shows the diagonal x matrix elements 
shared by a bare JPM and a JPM experiencing energy relax- 
ation. Energy relaxation matrix elements are represented by 
cirlces on the plots. 

As in the case of pure dephasing, it is instructive to see 
the effect of energy relaxation on the probability of pho- 
ton detection for specific input states. Figure 8 shows the 
detection probability as a function of time for both co- 
herent states and one photon Fock states. The short time 
oscillatory behaviour of the bare JPM detection proba- 
bility is strongly suppressed by energy relaxation because 
of the effective dephasing rate T 2 = 27^. In addition to 
this dephasing effect, the detection probability of a JPM 
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FIG. 7: Diagonal ^-matrix elements aT = Xj-lj+rj-ij+r 
corresponding to photon detection after loss of r photons. 
These elements are zero unless an energy dissipation mecha- 
nism is present. 

experiencing energy relaxation asymptotes more quickly 
to a smaller value than that of a bare JPM. 

A plot of the asymptotic value of detection probability 
as a function of energy relaxation rate, , is shown 
in figure 9 for a one photon Fock state input. The JPM 
experiences competing decay channels (energy relaxation 
and incoherent tunnelling into the measured state), only 
one of which results in photon detection. This reduces the 
time required for the probability to reach its asymptotic 
value as well as the overall probability that a detection 
will occur. In the case where these two decay rates are 
equal, the single photon detection probability will be ex- 
actly half of what it would be for a bare detector, as seen 
in figure 8. 

4- Dark Counts 

The final modification of interest is tunneling out of 
the |0) d state, a detection event which does not change 
the number or excitations in the cavity and is therefore 
considered a dark count. This is described by a master 
equation including the Lindblad operator Jo of equation 
(5). The x 1 -matrix has, in addition to all the non-zero 
elements of the bare JPM, additional non-zero elements 

4^0Vj,fc£{0,l,.JV-l}. (24) 

These elements correspond to detection events that occur 
without changing the fc-photon Fock state in the cavity. 
The JPM experiencing dark counts will appear to have a 
higher probability of photon detection than a bare JPM; 
however, this increased probability is due to false detec- 
tions. Dark counts limit the detector contrast by the 
ratio between true and false detections. 

We are not aware of any simple way to correct for all 
the effects of dark counts on detection probability with- 
out a priori information about the detected state; how- 
ever, by the appropriate choice of experimental parame- 
ters, the dark count rate can be made to be quite small 
for a JPM - as much as 2 to 3 orders of magnitude smaller 
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FIG. 8: Figure (a) shows the detection probability for a bare 
JPM and a JPM experiencing energy relaxation for a coherent 
state input (same as in 5(a)). Figure (b) shows the detection 
probability for a bare JPM and a JPM experiencing energy 
relaxation for a one photon Fock state input (same as in 5(b)). 
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FIG. 9: This figure shows the asymptotic value of detection 
probability for a JPM experiencing energy relaxation as a 
function of the energy relaxation rate T^ 1 . The black circles 
represent numerically simulated data points, while the red 
curve is a linear fit between the simulated data points. As 
expected, it is a monotonically decreasing function. 



FIG. 10: Detection probabilities for a bare JPM and one ex- 
periencing dark counts for (a) a one-photon Fock state input 
state and (b) an a = 1 coherent state input state. 



than the excited state tunneling rate [40]. It is therefore 
not unreasonable to simulate a dark count rate of 5 % 
of the excited state tunneling rate as a conservative esti- 
mate. 

Figure 10 shows the detection probability for a one- 
photon Fock state input with a 5% dark count rate, which 
is very similar to that of a bare JPM; however, for a 
coherent state with (n) = \a\ 2 = 1, dark counts signifi- 
cantly change the detection probability since the coherent 
state has a significant vacuum component. This devia- 
tion decreases as \a\ of the coherent state increases. 



D. Analytical Solutions in the Low T2 Regime 

While a full analytic solution of the system's master 
equation does not promise more illumination than the 
numerical results presented above, here we obtain the 
short T2 behavior of the detector by making appropriate 
approximations. We begin by defining the block cavity 
matrix 



pijit) 



(25) 



Id > \ / U/d 

where \i) d = I c (£) \i) d project £(t) on the basis states 
of the noninteracting JPM. We are interested in the un- 
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renormalized state of the cavity after photon detection, 
which corresponds to the cavity state p mTO (i). From the 
full master equation, it can be shown that 

Pmm{t) = (m| d e st £(0) |m) d = 7 ipn (t) (26) 

in the case that 70 = 0. Thus, the instantaneous time 
evolution of p mm depends only on pn, the unnormalized 
state conditioned on the JPM being in the excited state. 

The time evolution of p\\ is governed by a system of 
four first-order operator differential equations: 



P11 = 

Poo = 

P01 = 

Pro = 



ig (poi« - 

ig (pooa f 
ig (piio- 



- apoi) 
a f Pio) 

- a f Pn) 
apoo) - 



7iPn 



71 + K 2 



(27) 



-P01 



71 + K 2 



"PlO 



where we have introduced the pure dephasing rate n 2 = 
ijr for notational convenience. This system can be re- 
duced to a single fourth order operator differential equa- 
tion in terms of pn(t), as shown in equation (Al). 

Now for simplicity, we consider the occupation prob- 
ability of the cavity state \n) (n\ c given that a photon 
detection has occurred, 



Pn(t) 



(n\ c p mm (t') \n) c dt', 



which from (26) can be expressed as 

P„(i)=7i / x{t')dt', 
Jo 

where we have defined the matrix element 
x(t) = (n\ cPll (t) \n) c . 



(28) 



(29) 



(30) 



The full master equation reduces to the fourth order dif- 
ferential equation in x(t), given in (A2). 

Before attempting to solve the somewhat cumbersome 
equation (A2), we make the simplifying assumption that 
T 2 is the smallest time scale of the system's evolution, and 



keep only the highest order terms in 



^ < 1. 

K2 



This 



gives 



,(4) 



+ 



K 2 X 
1 



(3) 



1 



+ ^(2) + 



7 iK 2 2! (1) + «25 2 7i(! + n ) x 



0, (31) 



which can be solved with Laplace transforms. Defining 
X(s) = £[x(t)], equation (31) in the Laplace domain is 



X(s) = 2g 2 (n + l)x s 



«2 - 37i 



where x Q = x(0) = (n\ c pn(0) \n) c . 

While still very general and valid for all input states, 
equation (32) is still somewhat unwieldy. However, in 
the limit of short T 2 , we can assume that coherent oscil- 
lations between the cavity and JPM become incoherent 
tunneling. Defining T a .^-s.<5, 7 to be the tunneling rate 
from state |a) c (g) |/3) d to state \"f) c \ <S> |<5) d , we take 



r n ,0^n-l,l 
Fn,l— >n,m 



r n - 



l,l->n,0 



Ang 2 T 2 



7i 



(33) 



Here r n> o^. n -i,i is the incoherent tunnelling rate from 
the cavity into the JPM when n photons are present, and 
rn-i,i^n,o the rate for the inverse process. Both rates 
are broadened by short T 2 , which can be understood in 
terms of the Purcell effect. If we consider only Fock state 
inputs, the occupation probabilities 

P n j(t) = (n\ cPii (t)\n) c jG{0,l,m} (34) 

of the cavity being in the n-photon Fock state and the de- 
tector being in state |j) d obey the Pauli master equation. 
Using the rate in (33), this simplifies to 

Pnfi = n-f 2 (Pn-1,1 - Pn.a) 

Pn,l = (n+ 1)72 (Pn+1,Q - Pn,l) - JlPn,l 
Pn,m = llPn.l (35) 

where 72 = (2g) 2 T 2 = {2g) 2 / k 2 . For an n-photon Fock 
state input, -P n ,o(0) = 1 and the total number of excita- 
tions in the system is fixed to n, so at later times only 
1,1 and P n -i lU1 are nonzero. 
The equations (35) can be solved to find the detection 
probability as a function of time (the details of this are 
shown in Appendix A), giving 



Pnfi, Pi 



Pn,m(*) = 1 



71 72" 



+ 



s 



+ 



where 



^- 7 i - 2 72 n ± y / 7i 2 +47 2 2 n 2 ^ 



(36) 



(37) 



We can distinguish two regimes for this solution. In the 
tunneling-limited regime, <§; 7 2 n we find s + ~ — 7i/2 
and s_ ~ —2-f 2 n and have 



P„, m (t) = l-e-^/ 2 + 0( 7 i/72n) 



(38) 



In the opposite regime, photon capture is the slower, lim- 
iting process (7! ^> "f 2 n). In this regime, s + = —^ 2 n and 
S- = —71 — j 2 n and we find 



p 



1 - e-^ nt + 0(7 2 n/ 7 i) 



(39) 



s 4 + k 2 s 3 + -(k 2 s 2 + k 2 7i s ) + K 2g 2 Ji(n + 1 



(32)The differences between these regimes is evident in figure 
i 3(a); with added dephasing, 71 3> 72 n, and the detection 
probabilities P n>m — a n reach their asymptotic values 
more slowly and with n-dependence. Without the added 
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(b) Intermediate Ti Regime 

FIG. 11: These figures show the detection probability of a 
JPM experiencing pure dephasing in the low and interme- 
diate T2 regimes for one and two photon Fock state inputs. 
The detection probability obtained from the analytical so- 
lution described in this section is compared to a numerical 
simulation (via the Liouville supermatrix approach) of the 



detection probability, 
ure (b) i = 71. 



In figure (a) 



To 



IOOOO71 and in fig- 



V. COHERENT STATE TEST 

In this section, we propose a test to determine whether 
a given JPM's back action is closer to the lowering op- 
erator or the subtraction operator in order to correctly 
characterize the JPM. Additionally, such a test would al- 
low us to examine the effects of energy relaxation and 
pure dephasing on the time evolution of the back action. 



The Dependence of Detection Probability on 
Coherent State Power 



Consider a coherent state expressed in the Fock basis, 

I a) = e ~ J2^Lo '^H \ n ) ■ -^ ^ s straightforward to cal- 
culate the detection probabilities for both the lowering 
and subtraction operator back actions, as a function of 
a: 



Lowering Operator 
flow(a) = Tr [a \a) (a\ 

Subtraction Operator 



(40) 



-Psub(a) = Tr 



1-e 



"W =|al 2 - 



(a\ Bl 

JaL 4 , 
2 



0(\c 



(41) 



As equations (40) and (41) show, the difference in the 
detection probabilities occurs only in the nonlinear re- 
sponse, i.e., terms of the order of |a| 4 . This is consistent 
with the observation that a and B have identical ma- 
trix elements up to 1 photon but are different at higher 
photon numbers. By measuring the detection probability 
for coherent states of varying power and examining how 
this detection probability scales with the power of the co- 
herent state, it is possible to characterize the back action 
with respect to the lowering operator and the subtraction 
operator. 



dephasing, 71 <C 72^, the rate at which the detection 
probabilities reach their asymptotic value is determined 
by 7! , and independent of n. 

Equation (36) agrees to a high degree of accuracy with 
the numerical simulations for a JPM experiencing pure 
dephasing with a very short dephasing time, T2. As 
can be seen in figure 11(a) for one and two photon Fock 
states, this agreement occurs at all times. Interestingly, 
even at longer T2 the analytical solution is still a good 
approximation to the numerical simulation. By design, 
the analytical solution ignores coherent oscillations be- 
tween the cavity and the JPM, and so in the long T2 
regime the analytical solution does not display the oscil- 
latory behaviour of the numerical solution, but instead 
describes its average behavior. This can be seen in figure 
11(b), for one and two photon Fock state inputs. 



B. Bare JPM Coherent State Test 

First we apply this test to our simulations of a bare 
JPM, recalling that we expect B 1 (i) ~ a at short times 
(see equation (17)) and B 1 ^) ~ B m at long times (see 
discussion in sect. IV A). The proportionality constants 
affect the detection probability and not the structure of 
the back action, so we remove them by renormalizing (40) 
and (41) as follows: 



-Plow (a) = 
-Pub (a) = 



Pi ow (a)P data (a = q ) 
-Plow (a = an) 

P su b(Q!)P da ta(Q! = Op) 

P su b(a = an) 



(42) 
(43) 



where P da ta is the simulated detection probability, and 
an is the smallest value of a with simulated data. The 
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a-scaling of the JPM can be more directly compared to 
those of these rescaled lowering and subtraction opera- 
tors. Experimentally, this is important as it accounts for 
calibration uncertainties that may occur, such as from 
attenuation or imperfect impedance matching. 

Figures 12(a) and 12(d) show the detection probabil- 
ities of a bare JPM as a function of the power of the 
coherent input state a for a gt m value of 0.126 and 12.6 
respectively. As expected, the back action is very close 
to the lowering operator at short times - only deviating 
slightly at high powers - and at long times the back ac- 
tion is very close to the subtraction operator. At short 
times, the deviation from the lowering operator at high 
powers can be explained by equation (17), which only 
predicts the back action will be proportional to a for 
yfngt <C 1. Thus, at high powers (large \a\ = y/n), cor- 
rections to (17) will become important, as is this case for 
measurement times long compared to i cr jt = (gy/n) ■ 
Figures 12(b) and 12(c) show the back action at inter- 
mediate times. While 12(b) shows behavior intermediate 
between the lowering and subtraction operators, figure 
12(c) shows that a-scaling can actually fall below that 
of the subtraction operator. This effect is a result of the 
additional dephasing incurred in measurement. 



VI. OPTIMAL REGIME FOR A JPM 

Albeit based on a similar circuit, the bare operation 
conditions for the JPM are different from those of a phase 
qubit, where long T\ and T 2 are highly desirable. Op- 
erating a JPM at extremely long T 2 leads to the phe- 
nomenon of oscillating detection probability (figure 4), 
sub-subtraction back action (figure 12(c)) and additional 
dephasing of the cavity (figure 2(b)). This additional 
dephasing is undesirable as it destroys coherences in the 
original cavity state, irreversibly reducing its off-diagonal 
matrix elements, hence limiting the information available 
in a repeated measurement. 

It is thus advisable to operate the JPM in the short T 2 
regime. However, the effective photon-detector transfer 
rate, eq. (36), should be much shorter than the deco- 
herence rate of the cavity, and this places a lower bound 
on T 2 . On the other hand, T\ -processes always limit the 
measurement fidelity and should be avoided. One way to 
achieve the limit of long T\ with short T 2 is to damp the 
JPM with a frequency-dependent impedance with low- 
pass character, e.g., along the lines of [71], by shunting 
the JPM with an ii?-element. 



VII. CONCLUSIONS 



C. The Effects of Pure Dephasing and Energy 
Relaxation 

We now study the affects of energy relaxation and pure 
dephasing on the a-scaling of the coherent state test. Fig- 
ure 13 shows the detection probability of a JPM experi- 
encing energy relaxation at a rate jr — 71 as a function 
of a, at the same times as those of figure 12. As can 
be seen in figure 13(d) (which represents the long time 
steady state of the back action), the major effect of en- 
ergy relaxation is to prevent the back action from fully 
transitioning to the subtraction operator at long times, 
but rather it asymptotes to an operator in the intermedi- 
ate regime. With the additional energy relaxation chan- 
nel present, the JPM detection probability becomes more 
sensitive to the number of photons in the cavity and does 
not fully approach the subtraction operator, which can- 
not resolve photon number. 

In addition, energy relaxation suppresses the sub- 
subtraction scaling at intermediate times (see figure 
13(c)). We expect this is due to the added dephasing 
at T2 = 2Ti, since figure 14 shows that dephasing alone 
suppresses the drop below that of the subtraction opera- 
tor at intermediate times (figure 14(c)). As can be seen, 
the effect of pure dephasing on a-scaling is similar to 
that of energy relaxation; however, instead of stopping 
the back action from transitioning to the subtraction op- 
erator, pure dephasing merely increases the time scale on 
which this transition occurs. 



In this paper we have analyzed the back action of a 
JPM on the microwave cavity state it measures. Numer- 
ical investigations of the cavity \ matrix conditioned on 
a detection event give us a convenient quantitative de- 
scription of the detection process while including several 
relevant environmental processes. At short times, the 
back action of a bare JPM is similar to the lowering op- 
erator, while at long times, its back action approaches 
the subtraction operator with additional cavity dephas- 
ing. This additional dephasing can be reduced by adding 
pure dephasing to the JPM, which dampens the coher- 
ent oscillations between the JPM and the cavity without 
compromising the purity of the cavity state. Energy re- 
laxation decreases the asymptotic value of the diagonal 
cavity \ matrix elements and the detection probability 
by a factor of 71/(71 + T\ X )- It is useful to develop a test 
to determine which regime the JPM is operating in for 
different measurement times, and the coherent state test 
is one such test that is straightforward to implement. 
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FIG. 12: For a bare JPM, the detection probability is shown for a measurement occurring at t m g values of 0.126 in (a), 1.26 in 
(b), 2.52 in (c) and 12.6 in (d). a runs from a = 0.01 to a = 1 in 0.01 intervals, and in each figure, all three curves are scaled to 
be equal at a = 0.01. The time in (d) is chosen such that it is representative of the long time steady state of the back action. 



Appendix A: Low T2 Analytics 



operator differential equation: 



The four first order operator differential equations of 
equation (28) can be reduced to the following fourth order 



+ (« 2 + 2 7l )^ + 



«2 + 71 W «2 + 57 M 2 

2D + 4g 



V +T) oV + 2^7, 



2 f k 2 + 7i 



(2) 



4 5 2 Po + 7i(^^)^o 



( K 2 + 27i)(2g 2 +2?o)+7i 
Pn — Ag^aa' puacr = 0. 



«2 +71 



o (1) 

Pll 



(Al) 



r 



The superoperator T>o[f] = g 2 {a)a, /} = g 2 {a*af+fa 'd) ourselves to the occupation probability of the cavity 
is introduced for notational convenience. We restrict 
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FIG. 13: For a JPM experiencing energy relaxation (with T\ — — as before), the detection probability is shown for a 
measurement occurring at t m g values of 0.126 in (a), 1.26 in (b), 2.52 in (c) and 12.6 in (d). a runs from a — 0.01 to a — 1 in 
0.01 intervals, and in each figure, all three curves are scaled to be equal at a = 0.01. 



state |n) (n\ (see equation (30)) and define x(t) 



(n\ c pn(t) \n) . Equation (Al) becomes 



i (4) + (K2+2rn)xW + 



(3) 



«2 + 7l \ f K 2 + 57 



4g 2 + 4. 9 V + 2 5 2 7l 



2 

«2 +7i 



+ 8g 4 n + 2g 2 njx 



+ 4g 2 n + Ag 2 x (2) 
K2 + 7i 



( K2 + 2 7l )(2 ff 2 + 2 3 2 n)+7 1 
-4g 4 (n + l) 2 a' = 



«2 + 71 



,(1) 



(A2) 



r 



This can be further simplified upon the assumption that 
Tjr is large, as shown in equation (31). 



To obtain equation (36), we rewrite the equations of 
(35) in matrix form 
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FIG. 14: For a JPM experiencing pure dephasing (with T2 = — as before), the detection probability is shown for a measurement 
occurring at t m g values of 0.126 in (a), 1.26 in (b), 2.52 in (c) and 12.6 in (d). a runs from a = 0.01 to a = 1 in 0.01 intervals, 
and in each figure, all three curves are scaled to be equal at a = 0.01. 



where the index indicating the photon number in the 
cavity has been suppressed. The Laplace transform of 
this system of equations is 



7 2 n 
-(71 + 72") 
7i 



where we have used the fact that Pi(0) = #io, and defined 
Vi{s) =C[P t (t)}. 

Solving the system of equations (A4) for "P m (s) gives 






V m {s) 



U7i72 



s (s 2 + s (2nj 2 + 7i) + ^7172) 



(A5) 



Using partial fractions and finding the residues of V m at 
the poles allows us to rewrite equation (A5) as follows: 



^ , n 1 7172" 

V m (s) = -+ n U 

s s+ — S- 



s + (s — s + ) s_(s— s_) 



. (A6) 

where s± are as defined in equation (37). The inverse 
Laplace transform of equation (A6) can easily be calcu- 
lated to obtain equation (36). 
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